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atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
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The  dimensional  stability  of  a  material  often  plays  an  important  role 
in  the  use  of  that  material  at  high  temperature,  particularly  in  structural 
applications.  The  determination  of  creep  behavior  at  high  temperature  Is 
frequently  complicated  by  a  nonuniform  temperature  distribution.  In  order 
to  approximate  the  ideal  situation  of  a  uniform  temperature  distribution, 
various  schemes  are  used,  such  as  reducing  the  cross  section  within  the  hot 
zone  to  restrict  the  highest  stress  to  a  suitably  small  region.  However, 
in  mauiy  practical  situations,  as  well  as  in  experiments  dealing  with  the 
measurement  of  creep  in  composites  (Ref.  1),  a  material  is  inevitably 
subjected  to  a  nonuniform  temperature  distribution.  The  total  creep 
elongation  is  then  a  sum  of  contributions  from  material  elements  with  a 
distribution  of  temperatures. 

In  this  report,  we  analyze  an  exaunple  of  creep  in  a  material  under  a 
nonuniform  temperature  distribution,  and  we  use  simple  but  realistic 
assumptions  for  that  distribution  and  for  the  material  creep  properties  as 
a  function  of  temperature.  This  analysis  was  motivated  by  measurements  of 
creep  at  high  temperatures  in  unidirectional,  fibrous  carbon-carbon 
composites  (Refs.  2  and  3) • 

A  sample  of  uniform  cross  section  was  used  because  of  ease  of 
fabrication.  The  sample  was  placed  in  a  furnace  with  a  central  peak 
temperature.  A  tensile  load  was  applied  to  the  sample  and  creep  elongation 
was  measured.  However,  creep  rates  versus  temperature  were  not  obtainable 
directly  from  the  creep  data,  due  to  the  temperature  distribution  across 
the  sample.  The  following  analysis  allowed  us  to  predict  both  creep  rates 
at  a  uniform  temperature  and  elongation  rates  under  other,  nonun  if orra 
temperature  distributions. 

The  sample  in  our  model  is  long  and  thin  (depicted  schematically  in 
Fig.  1),  having  length  2L  and  uniform  cross-sectional  area  A  in  the  high- 
temperature  region.  We  assume  that  all  strains  are  small,  e  <<  1,  and  that 
the  sample  is  under  tension,  P,  or  uniaxial  tensile  stress,  o  =  P/A.  Creep 
rate,  ^x),  cind  the  temperature  distribution,  T(x),  are  functions  of 
position  along  the  sample  axis,  x. 
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CREEP  RATE  PROFILE 


Fig.  1,  Schematic  of  idealized  experiment. 


The  average  creep  rate  over  the  sample  length  is  given  by 
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We  assume  the  material  has  a  simple  Dorn-Weertman  (Ref.  4)  temperature- 
dependent  creep  behavior: 


£  =  £^(0)  exp(-E/kT)  (2) 

where  EqCo)  is  a  rate  coefficient,  c  is  stress,  E  is  activation  energy, 
k  is  Boltzmann's  constant,  and  T  is  absolute  temperature.  A  simple  expres¬ 
sion  for  the  average  creep  rate  under  a  nonuniform  temperature  distribution 
is  obtained  by  inserting  Eq.  (2)  into  Eq.  (1): 
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We  also  assume  the  parabolic  temperature  distribution  (Refs.  3  and  5): 
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with  a  peak  temperature,  .  For  convenience,  the  temperature  at  the 
ends  of  the  sample  is  taken  to  be  absolute  zero,  which  is  a  reasonable 
approximation  because  the  contribution  to  the  integral  from  the  regions  of 
the  sample  near  the  ends  is  small.  We  define  a  =  E/kTjjj^jj  to  be  a  dimen¬ 
sionless  activation  energy,  and  eliminate  x  in  favor  of  q,  a  dimensionless 

1  /P 

variable,  in  Eqs.  (3)  and  (4)  through  the  substitution  x/L  =  [1  +  (a/q)]"  ^ 
Using  the  fact  that  T(x)  is  an  even  function  of  x,  we  then  obtain 
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This  integral  can  be  easily  evaluated  by  noting  the  results  of  experiments 
on  high-temperature  creep  in  carbon-carbon  composites  (Ref.  2).  Activation 
energies  of  the  order  of  100  kcal/mole  were  obtained  at  temperatures  of 
around  2500  K,  thus  giving  values  of  a  of  -20.  Consequently,  we  can  expand 
the  binomial  in  the  integrand  in  powers  of  q/o  and  integrate  term  by 
term.  Because  the  integral  is  over  all  q  >  0,  we  integrate  the  series 
outside  its  radius  of  convergence,  which  results  in  an  expression  for  the 
integral  that  is  an  asymptotic  (rather  than  convergent)  series  (Ref.  6)  in 
powers  of  1/a: 
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(6) 


An  alternative  approach  to  evaluating  the  integral  in  Eq.  (5)  is  the  method 
of  steepest  descents  (Ref.  6),  an  approach  that  gives  the  first  term  in  the 
series  in  Eq.  (6). 

For  experimentally  relevant  values  of  a,  the  series  in  Eq.  (6)  gives  a 
good  approximation  to  the  integral,  to  within  a  few  percent,  with  only  the 
first  two  terms.  For  a  value  of  a  -20,  it  can  be  seen  from  Eq.  (6) 
that  the  first-order  term  in  1/a  gives  a  correction  of  less  than  k%.  This 
suggests  in  fact  that,  to  a  reasonable  approximation, 
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Using  Eq.  (7)  to  analyze  the  elongation  rates  at  different  peak  tempera¬ 
tures,  we  can  obtain  the  activation  energy,  E,  and  pre-exponential 


In  general,  best  precision  is  obtained  with  the  series  by  locating  the 
minimum  term  and  truncating  the  series  one  term  before  it,  see  Ref.  6. 


The  result  in  Eq.  (7)  can  be  viewed  in  terms  of  an  "effective  gauge 
length,"  defined  by 


<e>L 


^eff  ■  ^^eff 


(8) 


where  e  is  the  strain  rate  for  a  sample  under  a  uniform  temperature  equal 
to  given  in  Eq.  (2). 

For  the  temperature  distribution  in  Eq.  (M), 
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Thus,  in  a  typical  experimental  example  (Refs.  2  and  3),  for  a  =  20  and  a 
sample  length,  2L,  of  91 .4  cm  (36  in. ) , 

Leff  =  ^  (^5.7  cm)  =  8.9  cm  (3.5  in.) 


Equation  (2)  can  be  used  with  to  obtain  e(2Lg^^),  which  is  the 

elongation  rate,  and  which  gives  total  elongation  when  multiplied  by  the 
time  interval. 

A  problem  mathematically  similar  to  that  of  calculating  creep  under  a 
nonuniform  temperature  distribution  is  determining  the  amount  of  diffusion 
or  degree  of  heat  treatment  when  the  sample  is  held  at  a  uniform  tempera¬ 
ture,  but  the  temperature  varies  as  a  function  of  time,  T(t).  For  example, 
following  Shewmon  (Ref.  7),  the  mean-squared  diffusion  distance  is 
expressed  as 


<x^>  =  2  D(T)t 


(10) 


where  the  diffusion  constant  D(T)  is  typically  of  the  form 


D  =  Dq  exp(-E/k:T) 
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If  temperature  varies  with  time  as  T(t),  then  the  total  mean-squared 
diffusion  distance  over  a  time  t'  is 

<x  >  =  2  J  D[T(t)]  dt  (11) 

0 

In  conclusion,  for  thermally  activated  creep  we  have  shown  that  an 
asymptotic  expansion  for  the  average  elongation  rate  can  be  used  to 
determine  the  creep  parameters  a  =  and  the  pre-exponential 

factor  £^(0).  We  have  done  this  for  a  parabolic  temperature  distribution, 
but  the  same  approach  can  be  used  for  other  temperature  distributions. 
Having  determined  the  parameters,  we  can  use  the  expression  in  Eq.  (3)  to 
predict  the  elongation  rate  for  any  given  temperature  distribution.  In 
addition,  we  have  indicated  that  analogous  problems,  such  as  heat  treat¬ 
ment,  are  of  similar  mathematical  form. 
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